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Abstract 

In this paper we discuss the dynamics and transport properties of a massive particle, in a time 
dependent periodic potential of the ratchet type, with a dissipative environment. The directional 
currents and characteristics of the motion are studied as the specific frictional coefficient varies, 
finding that the stationary regime is strongly dependent on this parameter. The maximal Lya- 
punov exponent and the current show large fluctuations and inversions, therefore for some range 
of the control parameter, this inertial ratchet could originate a mass separation device. Also an 
exploration of the effect of a random force on the system is performed. 
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I. INTRODUCTION 



Transport phenomena play a crucial role in many problems from physics, biology and 
social science. In particular, there have been an increasing interest in transport properties of 
the Brownian ratchets (BR), devices out of thermal equilibrium in which a nonzero net drift 
velocity may be obtained from fluctuations interacting with broken symmetry structures |l, 
Q]. These devices that where first proposed by Smoluchowsky [3^ and later discussed by 
Feynmann J] have deserved a great deal of attention in the literature (for a review see 
Ref. [5I). There is a wide diversity of areas in which BR's are applied, for instance to the 
study of molecular motors , the description of ion channels and molecular transport within 
cells Q, the treatment of Parrondo's paradoxical games j^, 9], optical ratchets and directed 
motion of laser cooled atoms jlo| and coupled Josephson junctions 11 1. In particular, these 



„o„e,u,.ibHum models become e^edaUy ...eresfug due to the. teehno.ogiea. aEpUcatio,. 



on nanoscales and microscales 
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13| such as microscopic particle separation |14j. llSl. Il6l|. 



In these models a net transport can be induced with fiuctuations associated with an 



additive force and noise, but also without noise, in overdamped [l7| or underdamped systems 
Q, [3]. In particular, inertial ratchets may exhibit a very complex and rich dynamics, both 
regular and chaotic motion. Moreover, they can display multiple reversal currents where 
their directions depend on the inertia term and the amount of friction. 

In this paper, the problem of transport in periodic asymmetric potential of the ratchet 
type is addresed. We study a model with an underdamped ratchet under the infiuence of 
spatio-temporal external perturbations without and with noise. Other authors, have studied 



18,hJ,M. One of 



the reversal current as a function of the intensity of an external driving 
the facts of this paper is to reveal the current reversal as a function of the specific frictional 
coefficient that provides a mass separator device. 

The present work is organized as follows: in Section [III we introduce the model under 
study; in Section UTTl the dynamical study of the system in particular in the transport regime 
is presented. Section HV] is devoted to the study of the role of noise to modify the char- 
acteristis of the dynamics and to induce the transport current . Finally in Section |V] we 
summarize our results. Some of technical details are given in an Appendix. 
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II. ANULAR INERTIAL RATCHET 



We consider non interactive massive (m) particles placed in a one- dimensional ring of 
radius R = 1, inmmersed in a dissipative environment, under a periodic potential and 
driven by an external force. The time evolution is given by the equation: 

mx = -fix - + F''{x,t). (1) 

ox 

In Eq. [T], < x{t) < 27r represents the coordinate of the particle, fi the damping coefficient 
and V{x)a is a one dimensional, periodic potential given by 



Kcos(^i^) , forO < X < 



-Vo COS [ \ ' -anj , for < x < a, 

fulfilling the periodicity condition Va{x + a) = Va{x) where a = 2tt/N being the number 
of wells (or sites) along the ring. The parameter a (a > 0) controls the left-right asymmetry 
of the potential. For a > 1(< 1) the minimum in each well of the ratchet is displaced 
towards the right (left) while a = 1 corresponds to left-right symmetric potential. Figure 
[D^a) displays this static potential for = 4 and a = 5. 

Particles are driven by the external periodic driving force F^^{x,t). We take this to be 
the gradient of a time dependent potential with a spatial periodicity that is twice the one 
of Va in order that consecutive wells alternate in time as absolute minima. This condition 
restricts N to be even. We therefore consider 

^Hrr ^ dV'^Ux,t) . , ,rfsin(A^x/2) 

F^^x.t) = -e = -esmiout) (3) 

that represents a longitudinal stationary wave along the ring. Time is measured in units 
of the period r = 2tt/uj of the external driving and e is the coupling strength. There is a 
critical value e = Ec such that if e < the whole potential (V^(x) + V'^^{x,t)) always has 
a minimum per site. For e > Ec the potential may loose minima corresponding to alternate 
sites as the time varies. Figure [Uj^b) displays the whole potential at r/4 for = 4 and a = 5 
and e < Ec- 

According the values of the parameters and the initial conditions the motion would be 
bounded (libration near the sites) or unbounded (rotations around the ring). 

To obtain a directional current in the unbounded dynamics, it is useful to study the 
symmetries of the equation of motion Eq.[T]. For hamiltonian systems (i.e. /x = ) with 
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FIG. 1: 

(a) Static potential with a = 5 as a function of the coordinate x. (b) Total potential 

Va{x) + V'^''{x,t) with e < ffc for t = 7r/2u; as a function of the coordinate x. 

a = 1, there are two symmetries that prevent directional transport: 

Si:x^x + {2k + l)'^,t^-t,v^-v; (4) 

271 

S2 : X {2k + l)— - X, t ^t, V ^ ~v. (5) 

5*1 symmetry leaves invariant Eq. [1] while S2 changes the sign the whole expression. In 
accord to both symmetries, for each solution with velocity v will be another solution with 
value —V. If a 7^ 1 S2, symmetry is removed but Si still holds. When yU > the only 
symmetry removed is Si. Therefore, for obtaining directional transport the values of a and 
H have to be a 7^ 1 and /i > 0. Taking this into account, we claim that the current is 
induced by damping. In the overdamped limit, such that the inertial term can be dropped 
out, the transport phenomena in the ratchet described by Eq. [1] (setting m = 0) has been 
well understood [21]. To determine inertial effects on the directional current, we will study 
the characteristics of the dynamics as a function of the mass. 
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For practical reasons, we have numerically solved Eq. [T] using the Fourier expansion of 
the potencial Va{x) Eq. [2] up to 20-th order rather than the original one (see appendix). 

III. TRANSPORT AND DYNAMICS. 

In this section we will discuss the characteristics of the directional transport and the 
dynamics of the ratchet as a function of the specific frictional coefficient fi/m. Let us 
consider the dimensionless version of Eq. [1] as: 



x = p 



-a; + 7 iF^iNx) - |- cos (A^x/2) sin (27rt) 



(6) 



where (3 = ^ is the specific frictional parameter, 7 = ^^^^ is the ratio between the 



external force and the dissipation, F{Nx) = and e' = e/Vo (in the following 

we drop the prime). 

All the calculations that we report where made for N = 4, Vo = 10, u = 6, fi = 1 (that 
is 7 = 407r/3) and e = 6.5 greater than Ec = 6.26007. We fix all the parameters except m, 
therefore we study the system as long as f3 is changing. 

We study the dynamical behavior of the system in the stationary regime starting from 
an ensamble of one hundred random initial conditions {Nin = 100) for each j3 . We have 
calculated the the mean velocity v, measure in site per period of the driving force, for each 
initial condition is: 

V27r/t-oo t ' ^ ^ 

where k is the winding number (that is the number of rotation that the particle makes 
around the ring). Then, the average ensamble velocity is given by: 

(^) = 7^E^.> (8) 

and its fiuctuations: 



, jsr. E(--. - (9) 



In the same way, the average maximum Lyapunov exponent is: 



{Lmax) = j;^^ L^max^ (10) 
-''in j—i 
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where L{^^^ corresponds to the maximum Lyapunov exponent for the j — th initial condition 



(calculated by algorithm given in Ref.|22|) and the fluctuation 



(Tit \ 



\ 



T.r y^(-^Tnaa: — {Lmax))'^ ■ (H) 
'in j—i 



In Figure [2] the magnitudes deflned above are displayed. By inspecting this flgure (a,b,c, 
and d), we found diverse and complementary information that can be extract from them. 
From Figure [2](a), we observe the average of the maximum Lyapunov exponent belonging to 
the mentioned stationary orbits. These exponents have positive and negative values. The 
negative values, as usual, would be related with regular behavior meanwhile the positive 
with irregular one, this assert holds when the dispersion shown in FiglJt^b) remains small. 
However from the same flgure, we found coexistence of different attractors. We discuss some 
of them, indicated by the arrows in Figj3l The average ensamble velocity is shown in Figj2](c) 
while in Figl2](d) its dispersion. It can be observe for /5 > 5 an average current {v)max = 2 
is established with zero dispersion. We remark that in such limit, th e sy stem is overdamped 



and the inertial term is negligeable for studying the current (see ref. 2l|]) althougth its effect 
is present in the Lyapunov exponent as can be seen from its fluctuations above f3 = 5 (FigjS]) 
related to the intrawell dynamics. For lower values of (3, a complex behavior is observed: 
the values of the current fluctuates from —2 to 2 in some regions while it vanishes in others. 
The fluctuations displayed in the flgures do not correspond to a statistical effect but to the 
complex structure of the dynamic (i.e. appearance or desappearance of diferent kind of 
orbits), what is therefore worth to a more detailed analysis. We varied (3 with a smaller 
step in order to get a zoom in a particular range as can be seen in FigjSl The arrows 
show particular values where there are coexistence of attractors. The label A correspond 
to a chaotic orbit (with almost no transport) together with a regular one whose current is 

V = —4/3. B is a region where arise a regular orbit of current v = 2 meanwhile the chaotic 
orbit persists near to (3 = 1.49. In C three regular orbits with current v = 2, v = —2 and 

V = are present. The label D corresponds to coexistence of two regular orbits of currents 

V = —2 and v = 2/3 respectively. 

Figure H] displays the position of the particle as a function of time and the projection of 
the trajectory onto the phase space. The cases shown are related to those cases discussed 
in the text. 
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FIG. 2: Average of the maximum Lyapunov exponent {Lmax) (&)> dispersion of the maximum 
Lyapunov exponent ai (b), average velocity {v) (c) and dispersion of average velocity cj^^,^ (d) as 
a function of the dissipation (3. The velocity is measured in terms of sites per period of the driven 
force. For sake of clarity in the four graphs the points are linked by lines. The vertical dashed 
lines limit the region shown in FiglSj 

To stress the coexistence of attractors for a given damping value, as discussed above, 
the histograms of the maximun Lyapunov exponent (^l,„„^) and the mean velocity {rijj) are 
displayed in FigJS] The highest value of ul^^^ in the row (c) is the sum of the values of riy 
for V = 2 and v = —2 because both orbits have the same maximal Lyapunov for this (3. 

It is interesting to note that the current (v) in Figl2t^c) shows abrupt inversions as a /? 
changes in a similar way to other ratchets when the parameters varies. As the variation 
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FIG. 3: Blow up of Fig. [2 

of P depends only on the mass m, similar particles but different masses will have different 
velocities. Therefore, this behavior originates a mass separation device. 

IV. THE EFFECT OF NOISE 



We now turn to study of transport current assisted by noise in presence of an exter- 
nal random force ^(t) fulfilling {C,(t)C,{t')) = a'^5{t — t'). We solve the following stochastic 
differential equation: 



x = (3 



-X + 7 I Fa(A^x) - - cos {Nx/2) sin (27rt) J + i{t) 



(12) 



In the noiseless overdamped limit (m — > 0), the motion is bounded {v = 0) or unbounded 
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FIG. 4: Projection of the trajectory in the phase space and the position of the particle versus time 
for differents values of /3. (a) Regular orbit with mean velocity v = 2. (b) Same as above with 
V = —2. (c) V = —4/3. (d) V = 2/3 and (e) chaotic behavior with v = 0. 

{v 7^ 0) according to e < Ec ot e > Sc- In both cases the stationary motion is a single 
regular orbit robust under shght variation of the coupling. As it was shown in ref. {21] the 
addition of an appropriate noise to the bounded regular orbit can induce a transition to the 
unbounded motion {v 7^ 0) rather than its destruction. In other words, the effect of noise 
on the bounded motion can mimic the increase of coupling e from a lower to a higher value 
than the critical, therefore we will have a brownian ratchet whose maximal current is near 
to those of the unbounded orbit. 

This scope would be quite different when the inertial term is non negligeable since the 
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FIG. 5: Histogram of the maximum Lyapunov exponent (left panel) and histogram of the mean 
velocity v (right panel) for different values of (3. All the calculations where obtained with an 
ensamble of 100 random initial conditions. Each row (from top to bottom) corresponds to the 
regions indicated by arrows A, B, C and D in Fig. [3l 

stationary unbounded motion is now more complex due to the fact that small variation in 
the coupling may introduce or destroy orbits (20']. On the other hand in the inertial system 
the bounded and the unbounded motion are separated by a fuzzy region of e where there is 
coexistence of both types of motion. 

Unlike the overdamped system the mere addition of noise to the bounded motion does 
not lead to obtain transport although it depends onto the robustness and the mean velocity 
of the unbounded orbits in the chosen parameter region (/3, 7, e) therefore it is not easy to 
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predict when there will be directional currents. As an example, in Figj6](a), we show the 
current obtained as a function of noise added to a bounded system {e = 4.1, (3 = 1.5708). In 
this case, the maximum value of the velocity corresponds to an optimal noise (near a = 3.5) 
but its value is far from the maximum value which is reached in the overdamped system. 

Now, we add the random force to a system where should have transport even in absence 
of noise (as that shown in [3]). The effect of noise will strongly depend upon the parameter 
(3 because, as we have already seen, it determines the properties of the transport in the 
stationary regime (current of the single orbits, coexistence of orbits and their stability etc.). 
Let us consider the deterministic system whose f3 = 1.987 is indicated by the arrow E in the 
Figl3J There are coexistence of two orbits. One of them, has current v = —2 and slightly 
negative maximal Lyapunov exponent and the other one has current v = —1 and slightly 
positive maximal Lyapunov exponent. In the first row of Fig. [7|the histograms of both are 
displayed. The mixing of two orbits ocurs in such a way that the total average current and 
the average maximal Lyapunov are (v) = 0.34, (Lmax) = 0.013 respectively. 

When the random force is added, for an optimal noise {a = 0.3), the stationary orbits 
have V = —2 and slightly negative maximal Lyapunov exponents. This can be observed in 
Figini^b) and (c) where {v) and {Lmax) as a function of noise a are displayed. For sake of 
clarity, in the second row of Fig|7|the histograms of v = —2 and {Lmax) (see Ref. |22|) are 
displayed when the noise is near of the optimal one. 



V. SUMMARY AND CONCLUSIONS 



In this work we have studied the dynamics in the stationary regime of an inertial ratchet 
as a function of the specific frictional coefficient {(3). This dynamics is strongly dependent 
on such coefficient since in some regions slight variation of this parameter can destroy sta- 
tionary orbits and/or it can create new ones. There are regions where two or three orbits 
of different characteristics coexist. This can be verify from the fluctuations of the maximal 
average Lyapunov exponent as long as the respectively histograms. This fact is reflected in 
the transport properties of the system, since the mean velocity of the orbits involved, should 
be quite different not only in their magnitude although in their directions. The large fluc- 
tuations and inversions showed in the average velocity when (3 varies are the manifestation 
of such phenomenum. This property could be useful to build a mass separation device. 
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FIG. 6: a) Average velocity {v) vs. a. The chosen parameters are £ = 4.1,a; = 6, /3 = 1.5708. The 
calculations where obtained with an ensamble of 200 random initial conditions for each value of a. 
b) Average velocity (v) vs. a. The chosen parameters are e = 6.5, oj = 6, P = 1.987. c) Average 
of the maximal Lyapunov (Lmax) as a function of a. The parameters are the same as in b). 

We have also explored the effect of noise in the system when a random force is added. 

In the overdamped limit, when the deterministic dynamics is robust enough against noise 
it is possible to obtain a directional current assisted by noise by adding a random force 
to the system when its motion is bounded. In other words, the presence of noise in mod- 
erate amounts induces the transitions bounded-unbounded (libration-rotation) before the 
destruction of the deterministic dynamics structure. If the inertial term becomes relevant, 
our results show that it is possible to obtain this effect but with a maximal current lesser 
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FIG. 7: The upper row displays the histograms of the mean velocity v and maximal Lyapunov 
exponent L^ax for deterministic systems for (5 = 1.987, e = 6.5. The lower row shows the same 
histograms when the random force a = 0.2 is added. 

than the ideal. This behavior may be due to two different causes: 1) The presence of the in- 
ertial term makes the deterministic dynamics less robust against the noise 2) The unbounded 
orbits where the system access by adding the noise have smaller values of v than the ideal 
V for the choosen parameters. This could be decided with a more exhaustive analysis which 
is beyond the scope of the present work. 

Next, we explore the inclusion of a random force in the deterministic system where the 
dynamics is already unbounded. We have consider a region where there are coexistence 
of two stationary orbits with different characteristics (in particular their v have opposite 
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sign such that the average {v) is shghtly positive). We observe, when noise is added, one 
orbit is not robust and it dessapears while the other one (with v = —2) persists. Then, the 
current of the system, not only has an inversion but increases its absolute value. We wish to 
emphasize that noise, in this case, stabilizes the dynamics as can be seen from the reduction 
of average maximal Lyapunov exponent (Lmax)- This is due to the fact that the destroyed 
orbit has positive Lyapunov exponent. 

Another possible effect of noise when it is added to the system in a region of (3 with 
coexistence of orbits, is to induce transitions between them (provided that both orbits are 
robust against noise). We have not observed this but a more comprehensive study should 
be make. 
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APPENDIX 



Taking in mind the future study of the quantum analogue of this kind of ratchet using the 
Bloch-Floquet formalism, we have considered the Fourier expansion of the Potential Eq. [2] 



rather than its exact expression. That is, for sites, a = and < x < 27r: 

Va{x) = VoJ2[M<^)cos{Nlx) + bi{a)sin{Nlx)] . 
1=1 

The coeficients ai{a) and bi{a) have annalytical expressions given by: 

2aZ'7r \ 



(A.l) 



ai[a) 



bi{a) = Fi{a 



Fi[a) sin 



1 + cos 



^1 + 
2aln \ 



where 



Fia) 



l + aj 
4/(a - l)(a + 1)3 



(A.2) 



(A.3) 



[(a + l)2-(2/)2][(2aO'-(a + l)']vr ' 
For (2/ — 1) = a or 1/(2/ — 1) = a these expessions are undeterminated but taking the 

appropiate limit, we obtain: 

, , 1 
21' 



ai[a 



k{a) = 0. 



(A.4) 



14 



As a mesure of the difference of expansion Eq. [A.l] up to 20 th order from the exact potential 
Eq. [2] we consider 

O = ^^^l^^(^-Pr-Ve.actydx = 6 X 10"^ (A.5) 

In Figure [T](a), the exact potential is undistinguishable from its Fourier expansion. 
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